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We construct all the bulk and boundary unitary cubic curvature parity invariant gravity 
theories in three dimensions in (anti)-de Sitter spaces. For bulk unitarity, our construction 
is based on the principle that the free theory of the cubic curvature theory reduces to one 
of the three known unitary theories which are the cosmological Einstein-Hilbert theory 
the quadratic theory of the scalar curvature or the new massive gravity (NMG). Bulk and 
boundary unitarity in NMG is in conflict; therefore, cubic theories that are unitary both in 
the bulk and on the boundary have free theories that reduce to the other two alternatives. 
We also study the unitarity of the Born-Infeld extensions of NMG to all orders in curvature. 
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I. INTRODUCTION 



In three dimensions, there are three parity invariant pure gravity theories that are known to be 
unitary in the sense of tachyon and ghost freedom at the tree level. These are the (cosmological) 
Einstein-Hilbert theory with no local degrees of freedom, the quadratic theory built from the 
curvature scalar with the Lagrangian density R — 2Aq + aR 2 which has a single massive scalar 
degree of freedom 38 1, and the new massive gravity (NMG) defined by the action [l|, 0] 



d 6 x 



aR - 2A m 2 + — ( W - 



(1) 
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that provides a nonlinear extension of the Pauli-Fierz massive spin-2 theory with two degrees of 
freedom. Here, a = ±1 or it could be set to zero to obtain a purely quadratic theory. The im- 
portant point is that, with some constraints on the parameters, these three theories exhaust the 
list of unitary pure gravity theories in (anti)-de Sitter [(A)dS] and flat spaces in three dimensions. 
Therefore, if one searches for a unitary theory built from arbitrary powers of the Ricci scalar and 
the tensor, then the propagator of that theory should reduce to one of these unitary theories, with 
possibly redefined parameters (such as mass, cosmological constant etc.). In flat backgrounds, the 
problem is trivial: Any higher derivative (cubic and more) deformation of the above theories is al- 
lowed since the propagators are intact in this background. But, in constant curvature backgrounds, 
which we shall deal with in this paper, generically, all the higher derivative terms contribute to 
the propagators and therefore the unitarity analysis is actually quite involved. However, as we 
shall show in detail, tree- level unitary theories can be constructed systematically by studying their 
propagators with the recently developed tools in 0] and with the earlier tools of [J] for analyzing 
the unitarity of a higher derivative theory around (A)dS backgrounds. In general, there are several 
motivations for introducing higher powers of curvature tensors in a gravity theory. First, string 
theory requires higher curvature corrections; for example, cubic curvature corrections are given 
in 0]. Secondly, in four dimensions, asymptotic safety approach to quantum gravity (see ja] for 
a review) involves contributions of the every possible term constructed by curvature tensors that 
is consistent with general covariance. Hence, in the effective field theory perspective, Einstein's 
gravity which is nonrenormalizable should be augmented with higher curvature terms obeying the 
symmetry of the theory. An efficient way of analyzing the effects of these higher curvature terms 
on the propagator structure, and consequently on the unitarity of the theory is considered in this 
paper. In fact, as an example, we will construct all the unitary theories in three dimensions that 
are built from at most the cubic powers of the Ricci tensor. 

Several extensions of NMG have already appeared recently: In 0], cubic and quartic extensions 
of NMG was found using the requirement that a simple (essentially integrable) holographic c- 
f unction exists. In @, 0], a Born-Infeld (BI) type action was defined which extends NMG up to any 
desired order in the curvature (and in particular reproduces the same cubic and quartic extensions 
of Fzf w ith fixed parameters at each order of the curvature) and which has a holographic c-function. 



In [10( , order by order extension of NMG was introduced again using the notion of a holographic 
c-function. This order by order extension also matches the curvature expansion of the Born-Infeld 
extended NMG @j. 

It is worth to stress again, in constructing a generic unitary theory at any powers of curvature, 
our main principle is the following: The propagator of the theory should reduce to the propagator 
of the known three unitary parity invariant theories after possible redefinitions of the parameters. 
Note that this principle is merely a restatement of the unitary extension of a theory and does 
not assume any strong conditions such as the existence of a simple holographic c-function or the 
condition that the resulting theory can be obtained from a Bl-type action. 

Up to now, we have discussed bulk unitarity only. For AdS spaces, unitarity on the boundary 
is also an important issue because of the AdS/CFT correspondence. Out of the three bulk unitary 
theories, NMG always gives a nonunitary theory on the boundary The other two theories have 
rather wide ranges of the parameters which allow both bulk and boundary unitarity. Therefore, in 
AdS, if a cubic theory is unitary in the bulk and on the boundary, then its free theory reduces to 
either cosmological Einstein-Hilbert or the R — 2Ao + aR 2 theory. 

The cubic theory found before 0, S] is a single member of the continuous family of bulk unitary 
theories that we shall present. Moreover, we will more directly show the region where this cubic 
theory is unitary. In principle, our analysis can be extended to any powers of curvature tensors 
and to any dimensions. We will also give two examples of arbitrary power theories: The so called 
Born-Infeld extension of new massive gravity and its close cousin [8[ , specifically we will show that 
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their propagators reduce to that of NMG. Namely, like the cubic theory found in 0], BINMG is 
unitary in the bulk only. 

Since NMG ([T|) plays an important role in the construction of cubic or higher order theories, 
let us recapitulate its properties. For proper ranges (which we shall discuss) of the dimensionless 
parameters a, Ao and the dimensionful parameter m 2 , NMG is a tree-level (bulk) unitary theory 
generically describing a massive spin-2 excitation with mass M 2 = (^— a + m 2 at the linearized 
level around both flat and (A)dS backgrounds fll. 13. [l~ll-[l5| . Here, the effective cosmological con- 
stant is A = Am 2 with A = —2 (<r ± y/1 + Ao). In what follows, we will work with the mostly plus 
signature, assume k 2 > 0, and our convention for the sign of the Riemann tensor follows from 
[Vp, Vj/] V a = R 17 V p . In flat backgrounds, unitarity analysis of this model is quite straightfor- 
ward and has been carried out in several places, but in (A)dS backgrounds the analysis is somewhat 
more complicated: In Q], the theory was shown to be formally equivalent to the Pauli-Fierz mas- 
sive gravity in (A)dS, and in [lj] direct gauge-invariant canonical analysis was carried out by 
decomposing the spin-2 field in its irreducible parts under the rotation group. 

The layout of the paper is as follows: In Section II, we start with the most general cubic action 
based on the Ricci tensor and the scalar, and find the equivalent quadratic action which has the 
same O (h 2 ) expansion, that is the expansion in metric perturbation, as the original cubic action. 
In Section III, we discussed the unitarity of Born-Infeld extensions of NMG. In the Appendix, we 
explicitly calculate the O (h 2 ) expansion of BINMG. 



II. UNITARY CUBIC THEORIES 



The most general cubic curvature theory built from the Ricci tensor and the scalar is 



I = — I d A x 



(2) 



where a, Ao, oj, 77, a, /3 and 7 are dimensionless parameters whose signs and numerical values are 
arbitrary at this stage except, we normalize a 2 = 1, and w 2 = 1 or w = 0. On the other hand, 
m 2 is of [Mass] 2 dimension and without loss of generality we choose m 2 > and k 2 > 0. In flat 
backgrounds, which necessarily requires Ao = 0, we know that for any a the theory is unitary only 
if lot] = 0. For u) = 0, the theory should have the "right" sign Einstein-Hilbert term with a = +1. 
Furthermore, if w is also set to zero in this case, then there is no propagating degree of freedom; 
while for 7 7 7^ there is a spin-0 excitation with mass m 2 = — > in order to have a nontachyonic 



behavior |13l . Il4j . For n = and u 7^ 0, NMG is recovered for a = — 1 with two spin-2 degrees 



of freedom having mass m 2 = with u > [l]. We will not consider the case when a = 0. 
Therefore, in flat space, the already known picture at the quadratic level does not change at the 
cubic or higher levels. Thus, the main question is to find possible ranges of these parameters for 
which this theory is unitary around its constant curvature vacua. To answer this question, one 
has to find the O (hft^j action where h^p = — g^ v and is the (A)dS vacuum (or vacua) 

for which = 2Xm 2 g IMU . One can directly compute the O {h 2 ^ action of but this is highly 
tedious and such a direct approach would be practically impossible for some arbitrary R n theories. 
Therefore, we will instead employ a technique developed in 0] which boils down to finding an 
equivalent quadratic action which has the same propagator and the same vacua. The procedure is 
quite effective and at no point one needs the complicated equations of motion. For more details 
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and uses of this technique see 0] • Let us now first find the maximally symmetric vacuum or vacua 
of ([2]). This can be done with the help of the equivalent quadratic action, as we just said, but in a 
simpler way the vacuum can also be found from an equivalent linear theory. This follows from 



d 6 x£(R,R lll/ ) = / d i xC[R,R lw )+ / d A x 



5C 



5g 



(3) 



9y.v 



+ 



d 3 x8g 



arf 



5C 



5g a P5g^ u 



s<r+ 



9ij,v 



where C = y/—gf (R, R^v), and by equivalent linear action we mean an action which has the 
same O (hP) and O (h) expansions as (J3J), and equivalent quadratic action has the same O(hP), 
O (h) and O (h 2 ) expansions as given in Q. To find the equivalent linear or quadratic actions, 
f{R,R l j iV ) should be expanded to linear (or quadratic) order in the curvature around (R, R^v). 
The important point is that from the linear (or quadratic) expansion in curvature one gets all 
the 0(h) [or O (h 2 ) ] terms oi f{R,R lw ). Therefore, the expansion in small curvature is not an 
approximation as far as the vacuum and the propagator of the full theory is considered. [In these 

expansions one has to keep in mind that O (h n ) terms come from the J2i=o (R ~ ^) expansions.] 
We can now start our computation and find the vacua of ([2]). One further simplification is to 
consider the Lagrangian density as a function of i?^, in order not to introduce the metric or its 
inverse during the expansion. Therefore, we have 



-R 2 ) + 



+ 



a 
6m 4 



(4) 



Then, expanding / (R£) to the first order around the yet to be found background ( R£ = 2Xm 2 5^ 
with the assumption of small fluctuations [that is ( Rg — Ra ) being small] as 



f(R») = f (Rt + 



Of 



dR c , 



R<$-R a p)+0 



T>a net 
Kp - Rp 



(5) 



one obtains the equivalent linear Lagrangian density <?ii n -equai (Rtf) after dropping the quadratic 
order as 



fflin-equal 



o\2 

-2A + — (w - 37?) - 8aA 3 (1 + 3/3 + 9 7 ) 



m 



a--(u-3rj) + 2aX 2 (1 + 3/3 + 97) 



R. 



(6) 



Therefore, the equivalent linear action becomes 



-^lin-equal 



1 



d 3 x^/^g 



a-^(u-3ri)+ 2a\ 2 (1 + 3/3 + 9 7 ) 



R _ [4A - 3 (a; - 3??) A 2 + 16«A 3 (1 + 3/3 + 9 7 )] 2 
[2ct- A(w-3r/) + 4aA 2 (l + 3/3 + 97)] m 



(7) 



Let us stress again that ([7J and ^ have the same O (h°) and O (h) expansions. Since O (h) expan- 
sion of (J?]) evaluated at g^ u just gives the equations of motion, that is the Einstein tensor evaluated 
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in the vacuum in this case, we can easily read the vacuum, by comparing it to y/—g (R — 2Xm 2 ) 
and find 



2A 



4A - 3 (w - 3??) A 2 + 16a A 3 (1 + 3/3 + 97) 
2a - A (w - 3t/) + 4aA 2 (1 + 3/3 + 97) 



4aA + A 2 (w - 3/7) - 8aA 3 (1 + 3/3 + 97) = 4A , 

(8) 



which has always at least one real root for generic values of the parameters: Therefore, unlike the 
NMG case which requires Ao > — 1 for (A)dS to be the vacuum, for any Ao , © nas a maximally 
symmetric vacuum. At this stage, no restriction exists on the ranges of the parameters, but as we 
will see now, unitarity of the theory will constrain some of these parameters. 

Let us now find the equivalent quadratic action by expanding / (R£) up to second order in the 
curvature: 



Squad-equal (K) = f [K) + 

where 

f(K) 
Of 



df_ 

dm 



Rr — Rr) + x 



d 2 f 



dR p a 8R c r _ 



Kp - Up 



R P a-K), (9) 



3A 2 



6aX - 2A - (u - 3r/) + 4aA 3 (1 + 3/3 + 97) 



m 



dRI 



a - - (w - 3??) + 2qA 2 (1 + 3/3 + 9 7 ) 



5? 

u r\ 1 



(10) 



d 2 f 
dR p adR c r 



pirn 



in. 



[u + a\(l + P)]6$SZ 



U -3" 



(2/3 + 9 7 ) 



Then, collecting all these we get the equivalent quadratic Lagrangian density 



Squad-equal (i#) = "2A + 4«A 3 (1 + 3/3 + 9 7 ) 



+ ^[w + aA(l + P)]R 



2 



m 2 + 



8m 2 



cr-2aA 2 (l + 3^ + 97) 



1 

UJ 77 

3 ' 



~9~ 



(2/3 + 9 7 ) 



R 



R- 



(11) 



whose O (h 2 ), O (h) and O (h°) expansions match the same expansions of ([2]). At this stage, it is 
clear that there are three different ways for the general cubic theory ([2]) to be unitary: Its equivalent 
quadratic action (jlip can be, with redefined parameters, equal to the cosmological Einstein-Hilbert 
theory or R + aR 2 theory or NMG. [Again, we exclude the case for which Einstein-Hilbert term 
drops out.] First, it pays to rewrite the equivalent quadratic action as 



quad-equal 



d 3 x 



aR - 2A 7n + 



777/ 



R 



7 ' v 



R 2 + 



8m z 



;R 2 



(12) 



where 



a = a - 2a\ 2 (1 + 3/3 + 97) 



A - 2aA 3 (1 + 3/3 + 97) 



cj = u + aX(l + P), 77 = 77 + — (9 + 25/3 + 72 7 ). (13) 

Here, it is worth restating that A appearing in the redefined parameters is the vacuum of 
satisfying ([5]) which can also be directly obtained by computing the vacuum of (|12p which reads 
from the somewhat simpler looking expression 

crA + l - (a) - 37]) A 2 = Ao- (14) 



Canonical analysis of (|12[) have shown that there are generically three, not necessarily unitary, 
degrees of freedom with the masses 1_4{ : 



m 2 



fj 2 V 3f)J 
Co 2 2 Co 



m 2 helicity-0 mode, (15) 
m 2 helicity- ± 2 modes. (16) 



For (j!2p to be unitary, the necessary but not sufficient condition is Cofj = which again exhausts 
all three unitary theories. Among these theories, NMG, for which fj = 0, seems to be the most 
interesting one with spin-2 excitations (scalar mode decouples), therefore we start with it. But, 
NMG in (A)dS is not unitary by default: There are constraints on the parameters which we discuss 
below. Since the parameters appear in certain combinations let us define £ = 2a (1 + 3/3 + 97) and 
X = ot (1 + j3), then the effective parameters (fTB"|) become 

a = a - A 2 £, A = A - A 3 £, 

Co = oo + X X , V = V + ^( X + 40- (17) 

A. Reducing the cubic theory to NMG in (A)dS 

Setting fj = 0, the equivalent quadratic action (|12|) reduces to NMG with m 2 = ^— £ + ^Xj m 2 

where A = — % i^a ± \j 'a 2 + CoX^j which requires a 2 + CoXq > 0. The theory is unitary if 

(Xm 2 - 2a^pj > which comes from the ghost freedom requirement of 0] and reduces to 
CoX — 2a > in our notation. This requirement can be seen by rewriting NMG in the form of 
a massive Pauli-Fierz theory at the linearized level. In de Sitter case (A > 0), there is also the 
Higuchi bound [3j m 2 > Xm 2 which becomes A < 0, and in anti-de Sitter case (A < 0), 

there is the Breitenlohner-Freedman (BF) bound [17( m 2 > Xm 2 which is exactly like the Higuchi 
bound for this three-dimensional case. [Strictly speaking BF bound was derived for massive scalar 
field in AdS, but it works for massive spin-2 field as well [3]] In this setting, unitarity analysis of 
(j!2H for fj = is the same as NMG with an essential difference: a and Co are not in general ±1. 
However, as implied by the unitarity constraints, unitary regions can be classified according to the 
signs of a and Co just like in the case of NMG. Since the unitarity regions of NMG in (A)dS were 
studied in detail in we will not repeat the analysis here, but simply give an example in AdS 
(A < 0). Choose a < and Co > 0: BF bound is automatically satisfied, so the unique constraint 

on the vacuum of the theory is A > ^ with A = — % i^a + \J a 2 + CoXq^ which can be achieved if 

the parameters of the theory satisfy the inequality 

< A < (18) 
00 

This is a rather weak condition on the parameters, therefore there is a continuum of unitary 
theories. 

1. Choose a = — 1 and 00 = 1: For the sake of simplicity, let us further assume 77 = which 
fixes £ = — 4 that yields 7 = — 25 y 2 t " 9 m terms of the original parameters of the theory (we 
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discuss j] 7^ cases below). Then, for Ao < there is no unitary theory, but for Ao > the 
theory is unitary if the following conditions are met: 



X < -, and < A < — 2 



I Y> -, and < Ao < — ) . 
4 X/ 

(19) 

For example, consider the x = case, it is unitary for < Ao < 3 with the same vacuum 
as NMG, A = 2 (1 — y/1 + Ao)- In fact, NMG with a = is a member of this family, since 
X = ol (1 + /?)• But, /3 = —1 gives a cubic order extension which is probably the simplest 
unitary one parameter extension of NMG with the action 



K Z 



d 3 x 



(20) 



+ 



a 
6m 4 



T>liv pap p pi i pc 

Or K u tl all — JrCK^ + -it 



with an arbitrary a. The other one parameter extension of NMG introduced in [3j is also a 
member of r) = and rj = family of unitary theories, for this case one chooses j3 = 
which then fixes 7 = 17/64 yielding an action 



d 3 x 



-R - 2Xr ) m 1 + 



R 



flu 



3 2 

77 R 



(21) 



4x 
3m 4 



Q 1 7 

Dfiu pap P P 2 i p3 



whose unitarity region is given in (|19p . [In fact, original sign choice for a is +1 in 0].] Note 
that for x = — 1/2, (12 ip reduces to the cubic order expansion of BINMG which is unitary 
for < A < -8 + 6y/3. 

Let us also give an example for rj 7^ 0. For simplicity choose £ = which yields i] = 
then choosing Aq = 1 yields the unitarity region — 3 < x < 1 f° r the theory 



I 



-R - 2m 2 + 



p2 



-R' 



Ax 
24m 2 



R\ 



(22) 



+ 



X 

6(1 + 0) 



{r? v R v a Rati + fiRRfu, (T^"^ 3 ) 



where /3 is arbitrary, and A is the vacuum of the theory. Let us stress that the propagator 
of this theory is exactly like NMG with redefined parameters. 



2. Choose a 
and 



T and 



T: Then, rj = theory is unitary if Aq > (Aq < is ruled out) 



X < and 



< Ao < 



1 + (1 - 4 X ) 



3/2 



6x 



X 



2X 2 



(23) 



For 77 7^ and with the choice £ = 0, the unitary region is Ao > and — 3 ^ < ^ ^ ~~ Jo' 

3. Choose a = 1 and u; = 1: Then, 7/ = theory has no unitary region. For ij 7^ 0, certain £ 
theories such as £ = 1 have unitary regions. 
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4. Choose a = 1 and oj = — 1: Then, 77 = theory is unitary if 

1 1 1 + 6y + (1 + 4y) 3/2 , N 

--<X<0 and -<A < ^— . 2 4 

For 77 7^ and with choice £ = 1, the unitary region is — 2 < % < for Ao = 1. 

The above discussion reveals just a sample unitary cubic theories. The other branches for various 
sign choices of a, oj, a, oj and existence or non-existence of 77 can be studied both in AdS and dS. 

Although classifying all the unitary theories of the form of ([2]) for all parameter choices is a 
tedious job, it is relatively easy to find the unitary regions if some parameters are fixed as in the 
cubic extension of NMG given in [3] and as in the case of BINMG 0, 0|. In existence of a 
holographic c-function in a specific form is the main theme, so in this AdS/CFT based context Ao 
is set to be negative Ao = — -k and c-function in the considered form can only exist, if j5 = — 9/8 
and 7 = 17/64 with an arbitrary a. Also, a = +1 is preferred, while oj is allowed to be both dbl. 
Then, the equivalent quadratic action becomes 



-^quad-equal — ^ \ A X 



(25) 

where the vacua of the theory satisfies 4 A + A 2 cj — f A 3 = 4Ao- The unitarity condition and the 
Higuchi/BF bounds in terms of the original parameters of the theory become Xoj — 2 — -j^p > 
and A + 2 ^8cj-qA) — ^' respectively- With this setting, the theory is unitary in AdS if a; = +1 
and a < -jp ^3Ao — 8 — (4 — Ao) 3//2 ) ; or if oj = —1, there are some constraints on a which are not 
particularly illuminating to write. For cubic order of BINMG, a is further set to be 4, but there is 
no unitary region for a = +1. On the other hand, for a = — 1 and oj = 1, cubic order of BINMG 
is unitary in dS if — 2 < Ao < and unitary in AdS if < Ao < 8 + 6\^3j. 

The above analysis shows that for nontrivial % (or a, (3 in terms of original parameters), there is 
generically a continuous family of unitary theories, and the cubic theory of [7H9J is just an example 
of this family. Just like in the NMG case, there are some special points which need further attention. 
For example, at m 2 = Am 2 a new scalar gauge invariance of the form 5f h !iv = Xm 2 g^ u (^ arises, 



and one has a partially massless theory with a single degree of freedom |19h21||. [Note that for the 
Pauli-Fierz spin-2 theory in (A)dS which is not a diffeomorphism invariant theory, at the partially 
massless point the new gauge invariance is of the form 5^h^ u = V^V^C + Am 2 ^^, but the higher 
derivative theories that we are dealing here are diffeomorphism invariant, and therefore, V^V^C 
part is simply part of the diffeomorphism invariance, and should not be counted as a new gauge 
symmetry] The theory defined by has unitary partially massless regions (in contrast to a point 
in NMG) for x°~ < J2 anci ^0 > — f (we have assumed rj = 0) with 



OJ 



A± = - (-2a ± ^4 + 3wA . (26) 



Another special point is A — 2^ = where m 2 = for which the linearized theory reduces to the 
Proca theory for massive spin-1 field which can be seen by first writing the equivalent quadratic 
action in the form of Pauli-Fierz action by use of an auxiliary field say /„ V) and then by integrating 

out the metric perturbation which then yields a massive spin-1 field with mass ^— 8?m 2 ^ . 

The details of this procedure has been given in [2]. An overall appears in the Lagrangian; 
therefore, for ghost freedom oj > 0, and hence a < is required for nontachyonic mass in the 
region ax > — \ and oj\q < 4 (we have assumed 77 = 0) with A+ = 4a + 2^4 — Aq for oj = 1 and 
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A_ = —4a + 2\/4 + Ao for u = — 1 in both dS and AdS. [NMG is unitary only in AdS for this spin-1 
limit.] 

In the above analysis, we required that the O (h 2 ) theory of ^ reduce to O (h 2 ) of NMG with 
redefined parameters. Next, we discuss the remaining two possibilities. 



B. Reducing the cubic theory to Einstein's theory in (A)dS 

Pure Einstein's theory in three dimensions is locally trivial. Namely, there is no propagating 
degree of freedom; but in any case it is a unitary theory, and therefore the cubic theory should be 
allowed to have the same O (h 2 ) form as Einstein's theory around (A)dS. This follows from (|12j) 
by setting the coefficients of R 2 and -R 2 „ to zero. One then obtains 

I = _L j d 3 x J=g~ {oR - 2A m 2 ) , (27) 

where 



a 



A 2 

Ao + -j (3r? - u) 



(28) 



with the vacuum A = 4 1 which reduces to A = <rAo (assume Ao / 0). Then, (3 and 7 can be 
determined in terms of other parameters in (|2|) as 



ft 



1 + 



UJ 



aa\o 



7 



2 3?? - 2buu 
9 ~ 72craAo ' 



(29) 



For unitarity, we should impose the right sign Einstein-Hilbert theory that is a 1 + ^ (3ry — uj) 



> 



0. Therefore, any cubic theory satisfying this constraint will be unitary, yet with no local degrees 
of freedom at the linearized level. As a simple example, consider oj = 0, r/ = 0, then one should 
have /3 = — 1 and 7 = 2/9, and a = +1 is required to have a unitary theory with the action 



1 



-9 



R - 2\ Q m 2 + 



a 
6m 4 



R^R, a R, 



Ctfl 



RR^ + 



-R 



(30) 



As in Sec ill A[ the cubic theory with arbitrary a and with choices (3 = — 1 and 7 = 2/9 turned out 
to be special. Actually, {R^R u a R afl - RRl„ + fi? 3 ) is the unique cubic curvature combination 
that does not effect the free theory in both fiat and (A)dS backgrounds. Let us give another 
interesting example in the case for cj 7^ for which the cubic theory 



= -2 / d 3 x ^g~\ oR - 2X m 2 + ( R^ v 



-R< 



(31) 



a 
6m 4 



ryfiU pap 



1 + 



aaXr 



25a; 

RR~,, + — I — 

^ V 9 72aaAr 



has the same O(h ), O(h) and O (h 2 ) expansions as (|27p . Although this theory involves two 
massive excitations in flat space; in (A)dS, there is no propagating degree of freedom. Unitary 

regions of (|3ip is given in Table I. In the r\ 7^ case, (5 and 7 are determined as f3 — — 1 ™ — ^ 

To have a unitary theory in AdS, rj < 



-1 'Y = -- - 

3^ constraint should be satisfied; while in dS one has 
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A 


<t(1-^) > 


w 


Unitary Region 


AdS 
o-Ao < 


<x = -1 


Ao >0 


uj\ > 4 


+1 


A > 4 


a = +l 


Ao <0 


cjA < 4 


-1 
+ 1 


-4 < A < 
A < 


dS 
o-Ao > 


<x = -1 


Ao <0 


cjAo > 4 


-1 


A < -4 


(T = +l 


A >0 


o;A < 4 


-1 
+ 1 


Ao > 
< Ao < 4 



Table I: Unitary regions for u> ^ and t] = 0. 



C. Reducing the cubic theory to R— 2Aq + aR 2 theory in (A)dS 

The third and the final option of how ^ can be unitary is that it has the same propagator as 
the R — 2Ao + aR 2 theory. For this to happen, the coefficient of R 2 in the equivalent quadratic 
Lagrangian density (fT2l) should be set to zero. Therefore, this determines [3 to be j3 — 1 qTT* 
Then, after using the vacuum equation 4<rA + A 2 (25w — 2>rj) + 8cvA 3 (2 — 97) = 4Ao, or in a slightly 
more efficient form 4<tA + A 2 (uj — 3rj) — 4£ A 3 = 4Ao, the equivalent quadratic action can be reduced 
to 



d s x y 



4a\ + A 2 (w - 3r]) - 8A 



m* + 4X °- X2 l"-^R. 
4A 



10 R2 < 32 > 



This theory is not unitary for generic values of the parameters. One-particle amplitude [13(] and 
the canonical analyses [l_4( of the action 



-(R- 2A ) + aR 2 

K 



(33) 



— where A is determined 



show that it describes a single massive excitation with mass m 2 s 

by A — Ao — 6aKA 2 = 0. For unitarity a > is required for both AdS and dS, and for dS mj > 
and for AdS we have the BF bound m 2 > A. Therefore, the mass of the scalar excitation described 
by (BSD is 



3A [12A - 8crA - A 2 (w - 3??)] 



16(cjA - A 



m 



0y 



(34) 



The analysis of the unitary regions follows similar to Sec JII Al above. We will not repeat the analysis 
in its full detail, but just give some examples of the regions where the cubic theory 



1 



d 3 x 



-g{ aR - 2X m 2 + 



R 



R 2 + 



8m 2 



R 2 



(35) 



a 
6m 4 



R^R^R, 



Ctfl 



1 + 



aaXo 



RR%+lI? 



that reduces to (|32|) is unitary or nonunitary. For concreteness, consider 7] = and uj = +1, then 
for a = —1, the theory is not unitary in dS. For a = +1, the theory is unitary if £ > j~ and 

In AdS, for a 



n3/2 



-1, the unitary region is £ < and Aq < For 



1 ^ \ ^ l+72g+(l+48g) J 
4? < A < 864|5 

cr = — 1, for any value of Ao there is a unitary region for £ < 0. The analysis for 77 7^ can also be 
done in the same lines. 
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D. Central charge and boundary unitarity 



In all the above analysis, we have considered bulk unitarity only. For the applications of 
AdS/CFT, boundary unitarity is also relevant. From the detailed work of 0], we know that for 
NMG bulk and boundary unitarity are in conflict. This conflict is not resolved in the cubic order 
extension Ml, or the infinite order extension of NMG The bulk and boundary unitarity 

conflict follows from the requirement that a positive central charge is not allowed for NMG in 
the region where NMG is bulk unitary. Therefore, it would be quite interesting to find both bulk 
and boundary unitary higher derivative theories. As we will see in this section, there are many 
such theories. First, recall that the central charge of a generic three-dimensional higher curvature 
gravity theory can be found by using [311434 1 



8tt 



\X\m 



dC 

9tlv dR^ 



(36) 



where the coefficient in front was put to conform to the normalization of Brown-Henneaux [35|] . It 
is easy to see that the central charge of a generic higher derivative theory can be computed directly 
from the equivalent quadratic action, since - is the first order term in the Taylor series 

expansion of the full Lagrangian around its constant curvature vacuum. This simple observation 
leads to a remarkable conclusion in the light of the discussion above: Any higher curvature theory 
that reduces to NMG cannot be unitary both in the bulk and on the boundary. This explains why 
an extension of NMG, be it cubic or any power, that has a free theory like NMG will not have 
unitarity on the boundary and in the bulk, and hence perhaps will not be relevant to AdS/CFT. 
But, any higher curvature theory that has the same free theory as the cosmological Einstein theory 
will be unitary both in the bulk and on the boundary. The theories constructed in Sec JIIBI have 
the central charge to be 

24 ^ (37) 



\/\X\K 2 m 

Both bulk and boundary unitarity requires a > 0. But, these are not the only theories that are 
unitary everywhere: Let us now consider the higher curvature theories that have the same free 
theory as the aR — 2Ao + aR 2 that we discussed in Sec JII CI The central charge of (flUj) with to = 
can be computed as 

24vr / 3nX\ , . 

A \K z m \ 2, J 



For unitarity fj > 0, and in AdS since A < we should have a > —^f- to have c > 0. We should 
check if this constraint is consistent with the other constraint (the BF bound) m 2 s > Am 2 with 
m 2 = (^jj — m? and the existence of a negative A satisfying the vacuum equation <tA— ^A 2 = Ao- 
One can find families of theories satisfying these bounds, let us give a simple example for which 
we take n = and oj = 1, then the action (|35|) is bulk and boundary unitary for 

, 24£ - (1 - 16f) 3/2 - 1 1 
a = +1 and £ < and ^ 2 < A < ^, 



.= -1 and --L<«<0 and _ ^ + (1 4- 1^ + 1 < ^ < _ 1 
lb 32^ z 4£ 

where £ was defined just before Sec ill Al 

To summarize, if a higher curvature theory is required to be unitary both in the bulk and on 
the boundary, then it should have the same free theory as either the cosmological Einstein-Hilbert 
theory, or the R — 2Aq + aR 2 theory with the constraints satisfying the bounds discussed above. 
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III. UNITARITY OF BINMG 



Up to now, we have constructed all the unitary cubic curvature theories in (A)dS. The procedure 
can be carried on to quartic or more powers of curvature, but here let us give two examples of 
Born-Infeld gravities which in principle include infinite powers of curvature. Our first example is 
the Born-Infeld extension of NMG was introduced in N| with the action 



-^BINMG 



a 



det [g + — 



det g 



(40) 



where G^ u = R^ u — ^g^R and a = ±1. This particular form of the action was chosen to reproduce 
the cosmological Einstein-Hilbert action at the first order in the curvature expansion and the NMG 
in the second order expansion. These two conditions are actually met by another Bl-type action 
that we shall discuss below which constitute our second example. On the other hand, the cubic and 
fourth order extensions of NMG given in [3| which was constructed with the help of a holographic 
c-function matches the same orders of (I4UI). C ertain aspects of BINMG such as its central charge 
jfj, Hf3|, c-functions 0], classical solutions |2d-[29| have been studied. We will study the unitarity 
of BINMG with two different methods: First, with the help of an equivalent quadratic action that 
we have employed above, and secondly we will explicitly calculate the second order expansion in 
metric perturbation hn V with the methods developed in [3J . These two methods obviously will give 
the same answer, but it is worth checking that the equivalent quadratic action method works with 
the help of the second more direct method for this infinite order theories. This more direct method 
is highly involved in terms of computation; therefore, we put it in the Appendix. 

Let us analyze the BINMG action by finding its equivalent quadratic action: To do that we 
have to expand the determinant in terms of traces which was done in [3o| 



-^BINMG 



4m 



d 6 x 



a 
2m 2 



a 



12m 4 



S 



1-* 

2 



(41) 



where K and S are defined as 



K = R 



-R 2 



S = 8R^ h R^ a R a v — QRR^ U + R . 



(42) 



The unique vacuum of (|4ip by directly studying the equations of motion was found in 0, as 



A = cxAn 1 - — 



An < 2. 



(43) 



In the spirit of the current work, let us verify this result by finding the equivalent linear action 
which circumvents the use of equations of motion. Let us define 



1 



2m 2 1 M v m 2 



1 / ^ 



12?n 4 



8R$RZRl - 6KR; (SIR*) + fan* 



1/2 



2 



(44) 



which assumes, as above, that R^ is the independent variable. Expanding / (R%) around its 
constant curvature background (r% = 2Xm 2 5^) to the first order in (r^ — R@\ as (0) one can 
find the equivalent linear Lagrangian density. For this one needs 



f{R<t) = (l-a\) 



3/2 



2 



(45) 
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which requires a\ < 1, 
df 



dm 



a 



P\m 4m 2 (l-o-A) 



3/2 



12m 4 



24R{r* - 12R.PR - 6RlR$6P + SR 2 l 



which requires crA ^ 1, then 



df 



ac 



4m 2 



%V1 - crA. 



(46) 



(47) 



With these results, the equivalent linear action for BINMG becomes 



-flin -equal 



cr 



VI -crX 



drx ^J^g < R - 4am z 



1 



1 + -A + 

2 2V1 -o-A 



(Ac 



where one can read the effective cosmological constant as 

1 



A = 2a 



1 + f A+ VT^I (A °- 2, 



2Vl - crA = 2 - A , 



(48) 



(49) 



which requires Ao < 2, and after taking the square of the equation, one obtains ( 143 p . 

Expansion of / (-R^) around the constant curvature background by using @ with the assumption 
of small fluctuations about the background requires the quantity 



d 2 f 



dRZdR* 



1 



2mVl - <rX 



5 l3 5 a 



> u a u p 



(50) 



Using this and (|45|). (|47j) : one obtains the equivalent quadratic action as 



<0(R 2 ) 



1 

7 2 



aR - 2m A X + — [Rf a ,-^R 



(51) 



where, for crA < 1, 



a 



a 



VI -<jX'' 



An — An — 2 + 



1 



1 



VI - a\ 



(52) 



Remarkably, the equivalent quadratic action turned out to be NMG with redefined parameters. 
Namely, the effect of all the terms beyond O (R 2 ) simply change the parameters of the O (R 2 ) 
expansion of the action which was NMG by construction. Let us stress again that this equivalent 
quadratic action has the same free theory (that is the propagator), same vacuum and same central 
charge as BINMG. Vacuum of BINMG in terms of the redefined parameters is 



CjX 2 + 4ctA - 4A 



0. 



(53) 



From the discussion in Sec JII A\ we know that NMG is unitary under two conditions CjX — 2a > 
and ^ + A < 0. Now, the question is whether these conditions are satisfied together with the 
BINMG condition Ao < 2 or not. A simple analysis shows that BINMG is unitary only for a = —1 
in AdS for < Ao < 2, and in dS for Ao < 0. Therefore, this analysis answers the question raised 
in 0] about the unitarity of the a = +1 theory in the negative. This is true for bulk unitarity, for 
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boundary unitarity recall the central charge from 0, [2^] , or just compute it from the equivalent 
action ([511 as 



31 
2G~z 



a 



00A 
~~2 



3a£ 
4G^ 



(2 - Ac) • 



(54) 



Since in AdS < Ao < 2, and a = —1, the theory is not unitary on the boundary just like NMG, 
or the cubic extension of NMG. The a = +1 theory is unitary on the boundary, but as we have 
just seen it is not unitary in the bulk. This is an expected result, because the free theory of 
BINMG is the same as the free theory of NMG with redefined parameters, and there is the obvious 
conflict between the bulk unitarity condition u\ — 2a > and the boundary unitarity condition 
2a - Cj\ > 0. 

We mentioned that there was a second Bl-type action that reproduces NMG in the curvature 
expansion. The action of this theory reads 



Am 



(Fx 



det 



<7 



R 



1 



Ao 
2 



V-detg , (55) 



which, by use of 



becomes 
4m 2 



drx 



det A = - 
6 



(TrAy - 3Tr^Tr (A 2 ) + 2Tr (A 



1 + 



2m 2 



R 



in. 



R 2 



1 



R 2 + 



3m 4 



/ 5 q 23 q 

[BTRfR^ - -RRl v + -R? 



(56) 



Ao 
2 



(57) 



Quite interestingly, this action reduces to NMG at O (h 2 ) with the same redefined parameters as 
the BINMG. Therefore, at the free level, these two theories cannot be distinguished. 



IV. CONCLUSION 

We have found all the unitary cubic curvature theories in three dimensions around constant 
curvature backgrounds. Without any further constraint, we have shown that unitarity in the bulk 
and on the boundary allows a large family of solutions as opposed to the cubic curvature theories 
that have appeared in the literature before, which allowed only bulk or boundary unitarity. The 
theories we have found should be studied in the context of AdS/CFT. We have also studied the 
unitarity of two Born-Infeld extensions of NMG which turned out to be unitary in the bulk only. 
Besides the parity violating extension with the addition of a Chern-Simons term and/or carrying 
out the unitarity analysis to O (R 4 ) , a quite physically relevant extension of our work is to find 
the unitary cubic curvature theories in four dimensions, which is currently under construction. 
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Appendix: O (h 2 ) Action of BINMG 



In this Appendix, we calculate explicitly O (h) and O (h 2 ) expansions of the BINMG action. 
First of all, let us find the constant curvature vacuum of (j4"U|) by explicitly calculating the first 
order action in the metric perturbation. In it was shown that O (h) of the generic Bl-type 
action 



d D x 



det (gu, v + Aa U ) - («A + 1) y/- det g 



(58) 



where A^ u is in the form A^ v = a [R^y + + O (R 2 ) with the definition = R^ — j^g^R 



is 



l O(h) 



where A$ is the first order term in the metric perturbation expansion of A^. Here, a is defined 
as A^ v = ag^v and for BINMG it becomes 



(1 + a) 



D-4 
2 



d D x ^ \(l + a) (g ptM A$ + h) - (1 + a)^ 1 (aA + 1) h 



(59) 



a 



~ ^9fiiyR ) = -oXg^ => a = -aX, 



(60) 



which, when inserted to the action, yields the constraint a > — 1 aX < 1. For BINMG, Aa V is 
Apv = ^ (RfM, - \gpvR) , then A$ and g pfJ/ A$ becomes 



where Bz v and i?/, are the linearized Ricci tensor and the linearized curvature scalar with the 
definitions 



1 /.=. 



Then, for BINMG with a = and k — t- k , the O (/t) action becomes 



(62) 



2m 2 



'O(h) 



K 2 y/TTa 
2m 2 



(fx x/-g (1 + a 



2m 2 



k 2 VTT> 



dTx sj-g 



1 + a) h - \/l + a 1 



R L + 2Xm 2 h 
X 



i - yj 

(63) 



/? 



then the constant curvature background equation of motion can be found as in (|43p from the 
coefficient of W . 

Now, let us turn to the explicit calculation of O (h 2 ) action for BINMG. In 0], the second order 
action in metric perturbation for (|58p in three dimensions was calculated as 



1 



l O(h 2 ) 



+h^{AW-^g^A$ 



kol\/\ + a 

1. 

'-jjLU 2 ' 



* v=i { ^ff^i) - \ (^4i ) ) 2 - ( x + fl ) ^45 

1- v / TT^(aA + l)] (/i 2 - 2/1^)1 . (64) 



16 



With the explicit form of A^ v for BINMG, let us calculate each term separately. First, the second 
line of the above equation takes the following form by use of the definition of the linearized Einstein 
tensor = R^ v — \g^ v RL — 2A/i /Jl/ in three dimensions and by use of the equation of motion; 



d 3 x 



1 - VI + a 1 



Ar 



a 

m 



Xvnf 



Secondly, let us calculate the terms quadratic in A™. There are two such terms AuoAf^ and 



• (65) 



Ig^Afiu J , and the first one becomes 



g^ - V^Vj, + 2Am g^j R L 



by using f d 3 x \/—gR\ and f d 3 x y '—gR^ ' R\ v which can be found as 

d 3 x y/^R 2 L = [ d 3 x V 3 ! \-h" v (g^U - V M V„ + 2Am 2 ^) R L 



(66) 



(67) 



/ 



- V M V^ + 2Am 2 g /iI/ ) i? L + (ug^ v - Xm 2 g llu R L 

(68) 



WXm 2 R^ + Xm 2 g fll/ R L + 12A 2 m 4 ^ 



where the background Bianchi identity and integration by parts have been used. The other term 
reads 



/ 



d 3 x V-i 



1 



16m 4 



4:Xm 2 g^ v R L - 4X 2 m A g IJiU h 



- (2) 

Let us consider g^A^J which is 



9 j vi/ o 



and using 



Am _ 
2 9 /Jos' 1 



and 



J d 3 x yTtsTR® = h» v \^piu + A™ V - ^-g^h 



(69) 



(70) 



(71) 



(72) 
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one gets 



J d 3 x yf-fjgr*$ = j d 3 x ypjjbr £^ + g^R L + Xm% v h - 2Xm 2 h 

This computation is somewhat lengthy, and one needs 

t°K% {BP vlia ) L = (~Rjl u + 3Am 2 V " Xm% v h 
and the two expressions involving linearized Christoffel connection whose definition is LT£ 



(73) 



(74) 



v^r^ (r^) L (ri) L = | d 3 x^ 



(75) 



(76)" 



Collecting all the terms and making use of the equations of motion, one obtains 



1 



2/eVl - <?\ 



d A x yf^gW 



(77) 



(a-3A)^ + ^ 



( 



g^U - V M V„ + 2Xm% u S j R L + (n^ - Xm 2 g^R L 



which can be compared to (25) of [23|]. Then, one can observe that this is the O (h?) of NMG with 
the redefined parameters given in 
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